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I present recent developments in the lattice QCD calculations of the light hadron spec- 
trum. Emphasis is placed on the limitation of the quenched approximation in reproducing 
the observed spectrum and indications that the discrepancy is reduced by introducing two 
flavors of light dynamical quarks. 

1. INTRODUCTION 

QCD is believed to be the fundamental theory of strong interactions. This recognition 
first came from perturbation theory, which enables us to describe hadronic processes with 
large momentum transfer. However, the conventional perturbative approach fails at low 
energies p<^ 1 GeV, where the strong coupling constant is of order unity. 

In this region, lattice QCD provides us with a non-perturbative tool to calculate physical 
quantities of hadrons from first principles. A lot of works have been devoted over two 
decades to numerical simulations of lattice QCD, and gradually revealed that numerical 
method is a powerful and practical tool. 

Lattice QCD has been applied to calculations of the hadronic spectra including glueballs 
and hybrids, the hadronic matrix elements of operators of importance to weak decays and 
other hadronic quantities such as the pion nucleon sigma term and the proton spin carried 
by the quarks. High temperature QCD has been an important research area of lattice 
QCD. The topological structure of the QCD vacuum has also been investigated, because 
it may provide a qualitative understanding of the mechanism responsible for confinement 
and spontaneous breakdown of chiral symmetry. 

Among various subjects pursued by lattice QCD, deriving the light hadron spectrum 
is a key step, because it will be a fundamental confirmation of QCD at low energies and 
give us confidence when we apply lattice QCD to calculations of other physical quantities. 
Considerable progress has been achieved in this field over the past two years 



2. SYSTEMATIC ERROR IN LATTICE CALCULATIONS 

The first attempt to calculate the light hadron spectrum was made in 1981 [0,0. Nu- 
merous studies devoted to this issue have revealed afterwards that we have to overcome 
several difficulties to obtain numerically precise results. 

Because of the limitation of computer power, most of works to calculate physical quan- 
tities of hadrons employ quenched approximation, in which we ignore sea quark effects. 
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Table 1 



Simulation parameters versus year. 


year 


1982 


1985 


1988 


1993 


1998 


Lt 


16 


24 


50 


50 


56-112 


Ls 


8 


~ 12 


16-24 


24-36 


32-64 


Lstt (fm) 


~ 0.8 


~ 1.2 


1.4-1.7 


2.0-2.5 


3.0 


a (fm) 


0.14-0.1 


0.14-0.1 


0.14-0.1 


0.1-0.07 


0.1-0.05 


rUg (MeV) 


> 120 


> 120 


> 40 


> 40 


160-30 


mps/my 


> 0.7 


> 0.7 


> 0.5 


^0.5 


0.75-0.4 


statistics 


20 


20 


50 


200 


150-800 



One may expect that the quenching error is not so large, because the valence quark model 
describes qualitatively the observed spectrum. However it it difficult to evaluate the mag- 
nitude of the error theoretically. In order to justify the approximation, we ffist have to 
estimate its magnitude by comparing the quenched spectrum with experiment. 

In lattice QCD, the prediction for the real world can be obtained only after we take the 
infinite volume and continuum limit. Therefore we have to control systematic errors due 
to finite lattice size and finite lattice spacing. In addition, because computational cost 
rapidly becomes very large as the quark mass decreases, one has to extrapolate results 
obtained at heavy quark masses to the physical up and down quark masses. Present 
typical calculations are still limited to the range of the quark mass 30 MeV. 

The ffist step toward the goal of the spectrum calculation is therefore to obtain the 
definitive result for the quenched QCD spectrum, by controlling all systematic errors 
arising from extrapolations mentioned above. 

3. QUENCHED QCD SPECTRUM 
3.1. A history 

Instead of reviewing a long history of quenched spectrum calculations, we reproduce 
in Table |^ typical simulation parameters versus year. Lt (Ls) in the table is the number 
of points in the temporal (spatial) direction and a is the lattice spacing. Before 1993, 
physical lattice sizes were limited to the range La^ 2 fm, which was smaller than a 
typical size of hadrons, and therefore works made before 1993 were in an exploratory 
stage. In the table, we reproduce, in addition to the quark mass, corresponding value of 
the pseudo-scalar vector mass ratio (mps/my), which is often quoted to show how close 
we are to the physical up and down quark masses, rrig ~ 5 MeV. The ratios were limited 
to the range mps/mv^O.5, {rriq'^ 40 MeV), which was far from mps/my = 0.176 at the 
physical point. 

A turning point of the calculation of the quenched spectrum was marked by the GFll 
collaboration around 1993 [Q. They made simulations at three lattice spacings ranging 
over a = 0.1-0.07 fm on lattices with La ^ 2.3 fm. After making continuum extrapolations 
and a finite size correction, they obtained results shown in Fig.0 In their calculation, 
masses of vr, p and K were taken as input to fix the lattice spacing, the physical (and 
degenerate) up and down quark masses and the strange quark mass. Although the entire 
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Figure 1. Quenched QCD spectrum re- Figure 2. The most recent result of the 
ported by the GFll collaboration [If. quenched QCD spectrum 



spectrum was not covered, their results turned out to be consistent with experiment, 
within one standard deviation, or 1-9 %. 

In spite of this encouraging result, there was no further work to give results for the 
entire spectrum with errors convincingly better than 5 %. In particular, the question of 
how the quenched QCD spectrum deviates from experiment was not answered. 

3.2. State-of-the-art calculation 

A report from the CP-PACS collaboration in 1998 changed the situation. They 
made an extensive simulation to obtain the quenched spectrum with errors of 1-2 % for 
mesons and 2-3 % for baryons. Their results have sufficient accuracy to discuss how the 
quenched spectrum deviates from the experimental spectrum. In Fig.^ are summarized 
their spectrum result. 

First, we notice that the discrepancy between the CP-PACS result and experiment is 
less than 11 percent. This means that the quenching error is of order 10 %, which is not 
so large as we have expected. This observation is important because we may expect that 
the quenching error is generally of a similar order also for other physical quantities. 

In Fig-H, one can see statistically significant and systematic deviation of the quenched 
spectrum from experiment, amounting to 7a for some particles. How the quenched spec- 
trum deviates is summarized as follows. If one uses K meson mass as input to fix the 
strange quark mass, 1) masses of vector mesons rriK* and are smaller than experiment, 
2) octet baryon masses are systematically smaller than experiment and 3) decuplet baryon 
mass splitting is smaller than experiment. If one uses m,^ instead of mx as input, mx* 
appears consistent with experiment and the discrepancies for baryon masses are much 
reduced. However, rriK turns out to be much higher. In other words, the meson hyperfine 
splitting remains smaller than experiment. 

In summary, the CP-PACS results show that the strange quark mass cannot be tuned 
in quenched QCD so that all strange hadron masses are in agreement with experiment. 
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Figure 3. Continuum extrapolations of me- Figure 4. Continuum extrapolations of 
son masses with used as input. Results baryon masses with itlk used as input, 
from two different chiral extrapolations are The CP-PACS results are compared with 



reproduced. See text for details. 



GFll's. 



3.3. Possible systematic error 

Controlling systematic errors is a basic element in spectrum calculations. We discuss 
in some detail the CP-PACS attempt to reduce them so that the reader can understand 
what was improved and how large the remaining errors are. 

They performed simulations on four lattices with La^ 3 fm. For these lattices, finite 
size effect is estimated to be smaller than 0.4 % for the worst case. 

The CP-PACS collaboration employed the Wilson's formulation of quarks on a lat- 
tice. For continuum extrapolations with the Wilson action, the leading correction to the 
continuum value is proportional to the lattice spacing: 

m(a) = m(a = 0) X [1 + Cia + (Csfl)^ H ]. (1) 

CP-PACS data are fitted well by a linear function of the lattice spacing, as shown in 
Figs. 1^ and ^. The coefficient of the linear term Ci ranges from to 0.28 GeV. One can 
evaluate the magnitude of higher order terms, assuming that Ci ~ C2. We find that the 
quadratic term has an effect of at most 1 % at a = 0.075 fm, the central value of the 
range of their lattice spacings. Therefore finite lattice spacing effect which may not be 
removed by linear extrapolations is estimated to be 1% or less. 

Chiral extrapolation is the most delicate issue in controlling systematic errors. The CP- 
PACS collaboration performed simulations at five quark masses corresponding to mps/my 
= 0.75, 0.7, 0.6, 0.5 and 0.4. Previous studies with the Wilson quark action were limited 
to the range mps/fny^ 0.5, or 40 MeV. Data at the smallest quark mass, rriq ?5i 30 
MeV, are indeed important for reliable chiral extrapolations. For example, the nucleon 
mass as a function of the quark mass shows a clear negative curvature, as shown in Fig.^, 
which may not be uncovered if the data at the smallest quark mass were lacking: Indeed, 
the GFll collaboration made a linear chiral extrapolation. Consequently, hadron masses 
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Figure 6. y = ^ps/i^qi + "^92) versus x = 
rUq^ + rUq^ . Filled (open) symbols represent 
degenerate (non-degenerate) data. 



at finite lattice spacing and hence continuum extrapolations are quite diff'erent from those 
of GFll. See Fig.^ for examples. 

Theoretically, chiral extrapolations in quenched QCD are more complicated than in full 
QCD. The quenched chiral perturbation theory (QxPT) ^ predicts characteristic singu- 
larities in hadron masses in the chiral limit. Therefore we first have to check numerically 
the validity of the prediction. The CP-PACS collaboration extensively investigated this 
issue and found several evidences for the existence of the quenched chiral singularities. We 
reproduce one example. In full QCD, the relation mpg oc holds. In quenched QCD, 
however, mpg has a logarithmic singular term and the ratio mp^/rriq diverges toward the 
chiral limit. As Fig.^ shows, the CP-PACS data exhibits a clear increase of the ratio. 

The CP-PACS collaboration employed QxPT mass formulae [|^^ for chiral extrapo- 
lations. They also repeated the whole analysis employing conventional polynomials in 
quark masses to investigate the effect of choosing totally different chiral ansatz. They 
found that the difference in the continuum limit was 1.5 % for the worst case, which was 
only a 1.5 (T effect. We reproduce their comparison of meson masses in Fig.^. 

From these considerations, we conclude that the CP-PACS collaboration successfully 
reduced systematic errors to the magnitude of their statistical ones. Because the deviation 
from experiment amounts to 4-7 a, taking account of systematic errors does not change 
their conclusions in any significant ways. 

Although the CP-PACS results are very convincing, a crosscheck is of course necessary. 
Two years ago, the MILC collaboration calculated the nucleon mass in the continuum 
limit using the Kogut-Susskind quark action, another popular formulation of quarks on 
a lattice. Their result rriN = 964(35) MeV is slightly larger than experiment, while 
rriN = 878(25) MeV from the CP-PACS is smaller. Understanding the difference, albeit 
only a 2.5 a effect, is important. Also, the MILC collaboration calculated only the nucleon 
mass. It is desired to establish the entire spectrum for the Kogut-Susskind action. 
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Figure 7. Comparison of rriN/mp at 
mps/mv =0.7 for the Wilson (W) and 
clover (C) actions in quenched QCD. 



Figure 8. mjv/mp in A^^ = 2 full QCD [|T|]. 
Abbreviations: P:standard, R:improved for 
gluons, W: Wilson, C: Clover for quarks. 



4. IMPROVEMENT PROGRAM 

The CP-PACS result above was obtained by a full use of one of the fastest computers in 
the world, CP-PACS, for about one year. Such a huge computational cost may make us 
pessimistic in doing realistic calculations in full QCD. Roughly speaking, computational 
cost for full QCD simulations is 100-1000 times larger than that for quenched QCD. Even 
in quenched QCD, there are many problems which need much more computer time than 
the simple spectrum calculation. 

A breakthrough to this situation may be given by improving lattice actions. The basic 
idea is as follows. Let us recall how we estimate the magnitude of the systematic error 
due to finite lattice spacing. Because the leading correction is 0(a) for the Wilson action, 
and we fit data with a linear function, the remaining error is of O(a^), which turned out 
to be about 1% at a ~ 0.075 fm. Now suppose that we can invent an improved lattice 
action which has in principle no 0(a) correction. If one uses the action and performs 
a fit using a quadratic function of the lattice spacing, the remaining error is reduced to 
0{a^). In this case, we can obtain values in the continuum limit with the same accuracy 
from simulations at a ~ 0.15 fm. This value is larger by a factor two than for the Wilson 
action. We perform simulations with the physical lattice size La being kept fixed, while 
computational cost depends on L, but not on La. Because computational cost becomes 
large faster than the number of lattice points L^, computational cost is reduced by at 
least a factor 16 for the improved action. 

Years ago, Symanzik argued a general method to improve lattice action, order by order 
in lattice spacing, by adding higher dimensional terms to the original one [|lOl- Sheik- 



holeslami and Wohlert then derived the 0(a) improved action [|TT| for the Wilson quark 
action. The form of the improved action is given by 

where g is the gauge coupling and F^j^ is the field strength. Because the lattice represen- 
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Figure 9. Continuum extrapolations of vec- 
tor meson masses in Nf = 2 QCD, com- 
pared to those in quenched QCD. 
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tation of F^^, looks like a clover with four leaves, the action is often called "clover action" . 
The coefficient 

Csw can be calculated perturbatively or non-perturbatively. 

Fig.|^ shows a compilation of the results for the nucleon to p mass ratio [mj^/mp) at 
mps/fny = 0.7 as a function of the lattice spacing. We clearly observe that improving 
actions significantly reduces scaling violation so that the ratio agrees with a phenomeno- 
logical value (the star in the figure) within 5% already at a ~ 0.4 fm, whereas the scaling 
violation for the Wilson action amounts to 20% at the same lattice spacing. 

Similar improvement is observed also for full QCD [|T2|. In Fig.|^ are compared rriN/mp 
at mps/mv = 0.8 in two fiavor QCD for four possible types of action combinations, the 
standard plaquette or an improved action for gluons, and the Wilson or the clover action 
for quarks. 

5. FULL QCD SPECTRUM 

The observation of improvement in full QCD encourages us to do a systematic study 
of the full QCD spectrum using improved actions. The CP-PACS collaboration is now 
working on this subject PJTB|, employing a renormalization group improved action for 
gluons and the clover action for quarks. They simulate two degenerate sea quarks, iden- 
tified with the up and down quarks, while the strange quark is treated in the quenched 
approximation. Calculations are made at three lattice spacings in the range a ~ 0.2 - 
0.1 fm, keeping the physical lattice size at La ~ 2.4 fm. For chiral extrapolations and 
interpolations to the physical strange quark mass, they use data at four sea quark masses 
in the range mps/my = 0.8-0.6 and five valence quark masses with m-ps/friv = 0.8-0.5. 

Our interest at this stage is how the discrepancy observed in quenched QCD is reduced, 
when one introduces two light dynamical quarks. We reproduce in Figs.^ and ^ the CP- 
PACS results as a function of the lattice spacing and compare them with the quenched 
spectrum. We clearly observe that vector meson masses extrapolate to values noticeably 
clover to experiment than those for quenched QCD. The remaining discrepancy might be 



8 



due to the quenching effect of the strange quark. Sea quark effects in baryon sector are 
less clear. We observe small scaling violation for baryons, which is an encouraging result. 
However, the mass results lie 5-10% larger than experiment. The discrepancy might be 
due to possible finite size effects {La ~ 2.4 fm for the full QCD data, while La ~ 3.0 fm for 
the quenched QCD data), or chiral extrapolation uncertainties (the smallest quark mass 
for full QCD corresponds to mps/my ~ 0.5, while ~ 0.4 for quenched QCD). Further 
study is necessary to clarify sea quark effects in baryon masses. 

6. CONCLUSIONS 

In summary, we now have precise results for the quenched QCD spectrum which re- 
produces the observed spectrum with discrepancies of order 10 %. Simulations with 
two dynamical quarks show indications that the discrepancies in quenched QCD for the 
strange hadron spectrum are reduced, although we still have to treat the strange quark 
in the quenched approximation. 

In my opinion, successful application of improved actions and fast development of com- 
puter power have combined to open a new era of lattice QCD simulations. A realistic 
calculation of full QCD will be achieved in the near future, which would shed new light 
on hadron physics at low energies. 

I thank members of CP-PACS group for their collaboration. I am grateful to Y.Iwasaki 
and A.Ukawa for valuable suggestions on the manuscript. This work is supported in part 
by the Grant-in- Aid No. 09304029 of Ministry of Education and University of Tsukuba 
Project Research. 
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